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Final  Report 

Proposal  Number:  54155-EV 

Quantifying  the  Stable  Boundary  Layer  Structure  and  Evolution  during  T-REX  2006 
Principle  Investigator:  Sen  Chiao,  Florida  Institute  of  Technology  and  San  Jose  State 
University 

1.  Statement  of  the  problem  studied 

A  comprehensive  study  integrating  surface  observations,  data  from  in-situ  measurements, 
and  a  nested  numerical  model  with  two  related  topics  was  proposed  in  this  project: 

(1)  Turbulent  Closure:  To  evaluate  the  two  primary  PBL  parameterizations  (Mellor- 
Yamada-Janjic  scheme  (MYJ)  and  Yonsei  University  scheme  (YSU))  as  well  as 
quantify  differences  at  a  fine  scale  model  output  using  the  different  turbulent 
mixing/diffusion  options  in  the  WRF-ARW  model;  and 

(2)  Terrain  Influences:  To  investigate  the  terrain  slope  effects  upon  the  horizontal  and 
vertical  eddy  viscosity  coefficients  calculation.  The  different  terrain  averaging  and 
filtering/smoothing  schemes  will  be  examined  in  order  to  quantify  their  impact  on 
SBL  forecast  results. 

The  working  hypotheses  were:  (1)  the  average  error  (bias)  of  modeled  sensible  heat  flux  at 
10  meters  is  within  10%  of  the  measured  flux  at  90%  confidence  level;  and  (2)  the  choice  of  sub¬ 
grid  terrain  representation  is  not  significant  in  the  error  of  averaged  winds  at  some  height  during 
stable  conditions. 

2.  Summary  of  the  most  important  results 

The  use  of  high-resolution  (i.e.,  grid  spacing  at  4  km  or  less)  modeling  is  essential  to  advance 
our  understanding  of  mountain-valley  flows  in  association  with  stable  boundary  layer  conditions. 
Observations  and  model  results  suggest  that  the  T-REX/EOP4  event  never  quite  achieves  a  truly 
steady-state  down-valley  flow  condition  near  the  surface,  as  fluctuations  in  the  jet  strength  and 
the  near  surface  mean  valley  wind  direction  appeared  to  occur  throughout  the  event.  More 
specifically,  surface  winds  with  a  distinct  and  persistent  westerly  component  were  seen 
throughout  the  night  and  early  morning  over  much  of  the  valley  floor.  The  simulation  results  also 
suggest  that  adding  two  additional  nests:  an  outer  13.5  km  and  an  inner  500  m  can  produce  a 
much  better  evolution  of  the  nocturnal  low-level  jet  and  especially  in  terms  of  the  speed  max.  It 
also  appears  as  if  the  addition  of  the  outer  13.5  km  nest  was  far  more  critical  then  the  addition  of 
the  inner  500  m  nest,  lending  support  to  the  idea  that  accurately  capturing  the  synoptic  condition 
was  a  critical  player  in  this  valley  low-level  jet  event.  Overall,  the  low-level  nocturnal  down- 
valley  jet  is  fairly  well  resolved  by  the  simulations  in  terms  of  timing,  spatial  location,  and  height 
above  the  valley  floor  (Chiao  and  Dumais  2013). 

In  T-REX  EOP4  simulations,  the  temperature  profiles  indicated  a  conditionally  stable 
environment  which  became  more  stable  as  the  surface  cold  pool  formed.  Estimates  of  the  depth 
of  the  cold  pool  are  about  100  m.  The  effects  of  this  on  mixing  within  the  boundary  layer  can  be 
seen  in  the  vertical  velocity  fields.  Stronger  vertical  motions  are  found  near  the  valley  slopes  and 
above  the  temperature  inversion.  The  motions  near  the  slopes  are  associated  with  the  nocturnal 
downslope  flows,  while  the  upper  level  motions  are  associated  with  detrainment  from  the  vertical 
wind  sheer  and  mixing  at  the  top  of  boundary  layer.  The  layered  flow  structure  also  follows  the 
vertical  velocities  pattern.  Slope  flows  dominate  the  surface  layer,  with  a  transition  to  valley  flow 
above  this  in  the  residual  layer.  Exchange  between  the  layers  is  implied  by  the  variation  of 
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vertical  velocities.  Moisture  appeared  to  be  entrained  from  the  surface  layer  to  the  residual  layer. 
The  externally  forced  flow  layer  (southerly  flow  near  3000  m)  also  shows  large  vertical 
velocities  in  the  boundary  with  the  valley  flow  but  these  do  not  persist  as  the  external  flow  grows 
in  strength.  Further  analysis  of  virtual  potential  temperature  profiles  as  well  as  flux  tower  data 
would  reveal  the  extent  of  mixing  and  exchange  between  layers. 

The  model  vertical  resolution  is  limited  by  the  complicity  of  terrain  within  Owens  Valley. 
The  lowest  vertical  level  is  at  a  height  of  102  m  above  the  ground,  which  meant  that  all  surface 
flows  are  not  explicitly  solved  but  resulted  from  similarity  theory  in  the  boundary  layer  scheme. 
Overall,  there  is  good  agreement  between  observations  and  the  model  for  stations  across  the 
valley  throughout  the  event.  Nevertheless,  there  is  no  relation  for  surface  air  temperature  in 
QNSE  scheme  similarity  theory  making  it  challenging  to  forecast  this  variable  given  model 
results  (Pattantyus  et  al.  2010a,  b;  Chiao  and  Dumais  2013). 

In  addition  to  the  stable  boundary  layer  research,  the  PI  and  a  graduate  student  (Travis 
Washington)  were  examining  the  sensitivity  of  planetary  boundary  layer  schemes  in  association 
with  landfalling  hurricanes  (e.g..  Tropical  Storm  Fay  2008).  The  goal  was  to  advance  our 
understanding  of  how  lower  boundary  forcing  affects  hurricane  intensity  before,  during  and  after 
landfall  (Washington  and  Chiao  2013).  Also,  the  PI  and  a  graduate  student  (Kim  Whitehall)  were 
investigating  the  island  scale  convective  initiation  and  its  sensitivity  with  boundary  layer 
parameterization  over  maritime  environment.  The  target  areas  include  Barbados  and  Puerto  Rico 
(Whitehall  et  al  2013). 

The  comparisons  of  the  WRF  model  planetary  boundary  layer  schemes  were  also 
conducted  to  study  a  downslope  windstorm  and  rotors  in  Las  Vegas  valley.  Two  events  (March 
20,  201 1  and  April  27,  2010)  occurred  near  the  Owens  Valley,  CA  have  been  investigated.  It  was 
one  of  the  objectives  that  the  high-resolution  model  results  can  be  evaluated  using  GOES-R 
baseline  products  such  as  the  Downslope  Mountain  Wave  Turbulence  (Reside  and  Chiao  2013). 
Using  Doppler  lidar  measurements,  a  cold  front  event  over  the  Bay  area  is  also  under 
investigation.  The  goal  is  to  better  evaluate  the  WRE  model  physics  (i.e.,  planetary  boundary 
layer  schemes)  under  urban  scale  condition.  Using  lidar  data,  a  detailed  boundary  layer  evolution 
was  able  to  record  (Elaiz  et  al.  2012).  In  addition  to  those  research  activities  listed  above,  two 
modified  PBL  schemes  based  on  YSU  and  MYJ  with  k-epsilon  closure  assumption  are 
developed.  Simulations  with  these  two  new  schemes  are  ongoing.  Results  will  be  compared  with 
data  collected  from  MATERHORN  field  experiment. 

Overall,  the  WRE  PBL  schemes  are  very  much  sensitive  to  complex  terrain  regions.  No 
single  parameterization  can  represent  all  events  (i.e.,  weather  conditions)  and  all  mountainous 
terrain  regions  (e.g.,  Sierra-Nevada).  The  comparisons  performed  for  terrain  center  around  a 
conventional  averaging  method  that  tends  to  reduce  mean  barrier  heights  and  fill-in  valleys,  and 
a  reflective  envelope  topography  method  that  essentially  results  in  much  steeper  model  terrain. 
Results  from  this  project  suggest  that  the  reflective  envelope  method  when  used  with  the  Quasi- 
Normal  Scale  Elimination  (QNSE)  and  Mellor  Yamada  Janjic  (MYJ)  schemes  reproduce  the 
surface  flow  the  best  in  the  Owen  Valley,  CA,  but  the  Yonsei  University  (YSU)  scheme 
performed  better  in  the  Las  Vegas  Valley.  The  inversion  height,  strength,  and  low-level  stability 
and  wind  speeds  near  the  surface  are  factors  to  be  further  investigated  (Pattantyus  et  al.  2011; 
Chiao  and  Dumais  2013).  Initial  and  time-dependent  lateral  boundary  conditions  may  be  further 
improved  via  data  assimilation  or  surface  obs.  nudging  (Eiserloh  and  Chiao,  2014). 
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